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ABSTRACT 

We investigate the modal properties of the r-modes of rotating neutron stars with the 
core filled with neutron and proton superfluids, taking account of entrainment effects between 
the superfluids. The stability of the r-modes against gravitational radiation reaction is also 
examined considering viscous dissipation due to shear and a damping mechanism called mutual 
friction between the superfluids in the core. We find the r-modes in the superfluid core are split 
into ordinary r-modes and superfluid r-modes, which we call, respectively, r°- and r®-modes. 

The two superfluids in the core flow together for the r°-modes, while they counter-move for 
the r'^-modes. For the r°-modes, the coefficient kq = limn^oi^/^^ is equal to 2m/[l'{l' + 1)], 
almost independent of the parameter rj that parameterizes the entrainment effects between the 
superfluids, where 11 is the angular frequency of rotation, to the oscillation frequency observed 
in the corotating frame of the star, and V and m are the indices of the spherical harmonic 
function representing the angular dependence of the r-modes. For the r'^-modes, on the other 
hand, kq is equal to 2m/[l'{V -|- 1)] at r; = 0 (no entrainment), and it almost linearly increases 
as 1 ] is increased from i] = 0. The r°-modes, for which w' = — v'^ oc correspond to the 

r-modes discussed by Lindblom & Mendell (2000), where v'^ and v'^ are the Eulerian velocity 
perturbations of the neutron and proton superfluids, respectively. The mutual friction in the 
superfluid core is found ineffective to stabilize the r-mode instability caused by the r°-mode 
except in a few narrow regions of rj. The r-mode instability caused by the r®-modes, on the 
other hand, is extremely weak and easily damped by dissipative processes in the star. 


Subject headings: instabilities — stars: neutron — stars: oscillations — stars : rotation 


1. Introduction 

One of the roles expected for the r-mode instability to play (Andersson 1998, Friedman & Morsink 
1998) is deceleration of the spin of newly born hot neutron stars by emitting gravitational waves that carry 
away the angular momentum of the star (e.g., Lindblom et al 1998). We know, however, that among older 
and colder neutron stars as found in LMXB systems there are many rapidly rotating neutron stars like a 
millisecond pulsar (see, e.g, Phinney & Kulkarni 1994). This fact suggests the possibility that the r-mode 
instability does not always work well to spin down the rapid rotation of the stars. For the r-modes in 
cold neutron stars with a solid crust, for example, Bildsten & Ushomirsky (2000) suggested a damping 
mechanism operating in the viscous boundary layer at the interface between the solid crust and the fluid 
core, to explain the clustering of spin frequencies around the value of 300Hz for accreting neutron stars 
in LMXB systems (van der Klis 2000; see also Andersson, Kokkotas, & Stergioulas 1999). For the modal 
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properties of the r-modes in neutron stars with a solid crust, see Yoshida & Lee (2001), who showed that 
the r-modes in the core are largely affected by resonance with the toroidal sound modes propagating in the 
solid crust. 

As neutron stars cool down below T ~ 10®K, neutrons and protons in the core are believed to be in 
superfluid states (e.g., Shapiro & Tuekolwsky 1983). In a rotating system of superfluids, it is well known 
that scattering between vortices in the superfluids and normal fluid particles produces dissipation called 
mutual friction (e.g., Khalatnikov 1965, Tilley & Tilley 1990). Therefore, for people who are interested 
in the r-modes instability, it was a serious concern whether the r-mode instability could survive the 
dissipation due to mutual friction in the superfluid core of cold neutron stars, e.g., in LMXBs. It was 
Lindblom & Mendell (2000) who first examined the damping effects of mutual friction in the core on the 
r-mode instability, and concluded that the mutual friction could not be strong enough to damp out the 
instability in most of the parameter domains which we are interested in. In their analysis of the r-modes 
in neutron stars, Lindblom & Mendell (2000) employed a perturbative method in which the spin angular 
frequency 12 is regarded as the infinitesimal parameter to expand the eigenfrequencies and eigenfunctions 
of the modes, and they looked for the r-modes with the scalings given by /?' = /r), — -I- me/ig/mp oc 

and w' = Vp — v'^ oc. where /ip, /i„, and /tg are the chemical potentials of the proton, neutron, and 
electron in the core, and Vp and are the velocities of the proton and neutron superfluids, and the prime 
(') indicates the Euler perturbation of the quantity. 

Recently, Andersson & Comer (2001) discussed the dynamics of superfluid neutron star cores, and 
confirmed an earlier result by Lee (1995) that there are no (/-modes propagating in the superfluid core. 
They also applied their argument to the r-modes in the core filled with neutron and proton superfluids, 
and suggested the existence of two distinct families of the r-modes in the core, i.e., r-modes for which the 
neutrons and the protons flow together, and those for which the neutrons and the protons are counter 
moving. Lindblom & Mendell (2000) considered the former family of the r-modes, which are less strongly 
affected by the mutual friction than the latter. 

In this paper, we employ a different method of calculation to investigate the r-modes in rotating 
neutron stars, although the basic equations describing the dynamics of superfluids in the core are essentially 
the same as those given in Lindblom & Mendell (1994). Our method of solution is a variant of that used in 
Lee & Saio (1986). Because the separation of variables is not possible for perturbations in rotating stars, we 
expand the perturbations in terms of spherical harmonic functions Y[^{9, cj)) with different Ts for a given m. 
We substitute the expansions into linearized basic equations to obtain a set of simultaneous linear ordinary 
differential equations of the expansion coefficients, which is to be solved as an eigenvalue problem of the 
oscillation frequency. In this method, we do not have to assume apriori a form of solutions for the r-modes 
in the lowest order of 12. In §2 we present the basic equations employed in this paper for the dynamics 
of superfluids in the core, and in §3 dissipation processes considered in this paper are described. §4 gives 
numerical results, and §5 and §6 are for discussions and conclusions. 


2. Oscillation Equations in the Superfiuid Core of Rotating Neutron Stars 

Microscopically, the superfluids in a rotating system move irrotationally everywhere except within 
the core of vortex lines. Averaging over many vortices in the fluids, we may define the average superfluid 
velocities < v >, which can satisfy the usual relation for uniform rotation Vx < v >= 20 in the equilibrium 
state (see, e.g., Feynman 1972). In the following we simply use v, instead of < >, to signify the superfluid 
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velocities. Hydrodynamic equations for a rotating superfluid based on the two-fluid model with the normal 
fluid and superfluid components are derived, for example, in Khalatnikov (1965). 

We derive basic equations governing superfluid motions in the neutron star core in the Newtonian 
dynamics, assuming uniform rotation of the star. The core is assumed to be filled with neutron and proton 
superfluids and a normal fluid of electron. We also assume perfect charge neutrality between the protons 
and electrons because the plasma frequency is much higher than the oscillation frequencies considered in 
this paper (see, e.g., Mendell 1991a). Since the transition temperatures Tc ~ 10®K to neutron and proton 
superfluids are much higher than the interior temperatures of old neutron stars (see, e.g., Epstein 1988), 
we assume that all the neutrons and protons in the core are in superfluid states and the normal fluid 
components of the fluids can be ignored. 

The basic hydrodynamic equations employed in this paper for the neutron and proton superfluids in 
rotating neutron stars are essentially the same as those given in Mendell (1991a) and Lindblom & Mendell 
(1994). In a fluid system in which two superfluids coexist, the entrainment between the two superfluid 
motions occurs because a Cooper pair of one fluid particles is affected by the force field produced by the 
other fluid particles (Andreev & Bashkin 1975 for a system of ^He and ^He superfluids). In the system 
of the neutron and proton superfluids in the core the entrainment effects between them are mediated by 
quantum mechanical nuclear force between the neutrons and the protons (e.g., Alpar et al 1984). Here, we 
introduce the entrainment effects in mass conservation equations, which are given by 

^ + V-j„ = 0, (1) 


and 




( 2 ) 


where pn and pp are the mass densities of the neutron and proton superfluids, and the mass current vectors 
and jp are defined as 

Jn — Pnn'^n Pnp'^pt ( 3 ) 


and 


jp 


PppVp + Ppn'^rij 


( 4 ) 


where and Vp denote the velocities of the neutron and the proton superfluids, respectively, and the 
coefficients p„„, pnp, Ppp, and Ppn are defined to satisfy /?„ = /?„„ -I- pnp and pp = Ppp + Ppn and pnp = Ppn 
under Galilean transformations. The mass conservation equation for the electron fluid is given by 


^ -b V • (peVe) = 0, 


( 5 ) 


where pe and Vg denote, respectively, the mass density and the velocity of the electron fluid. The velocity 
equation of the neutron superfluid is in an inertial frame given by 


-b Vn ■ VU„ = -V{pn + 5*) + — {Vp - V„) X (V X V„) 
Ot Pji 


( 6 ) 


where p„ is the chemical potential of neutron per unit mass, and 4^ is the gravitational potential. The term 
proportional to Pnp on the right hand side of equation (6) represents a drag force between neutrons and 
protons. If we assume perfect charge neutrality of the proton and electron plasma, we may have 


jplPp = '^e, and pplmp = pelme, 


( 7 ) 
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where rrip and rrie are the proton and the electron masses. The velocity equation for the proton-electron 
fluid is then given by 

771 \ 777 ( 777 \ O 

Vp H- -Ve + • VVp H- -Ve ■ VVg = -V /Xp H- -He + - — i^p “ fn) X (V X Vp) , (8) 

rup J Vfip \ rup J Pp 

where C = 1 + me Imp, and pp and pe are the chemical potentials per unit mass for the proton and electron, 
respectively. Note that we have neglected the entropy carried by the normal fluid of electron for simplicity. 
The Poisson equation is given by 

= 4ttGp, (9) 

where p = Pn + Pp + Pe and G is the gravitational constant. 

To linearise the hydrodynamic equations for the superfluids in rotating neutron stars, we assume that 
the neutron and proton superfluids and the electron normal fluid in an equilibrium state are in the same 
rotational motion with the angular velocity Q around the axis of rotation, which is along the z-axis. In a 
perturbed state, however, the neutron and the proton superfluids move differently from each other in the 
core, obeying their own governing equations. The mass current vectors are linearized to be 

Jn = P'nVo + j'n, and Jp = PpVQ + j^, (10) 

where the prime (') indicates the Eulerian perturbation of the quantity, and Vno = VpQ = Vq = rsin^fle^ is 
the fluid velocity in the equilibrium state, and 

jn ~ PnnVj, PnpVp, and = Pp^V^ -\- PppV^ (H) 



are the perturbed mass current vectors in a corotating frame. The perturbed superfluid velocities are then 
given in terms of and jp as 


where 


and 


= M;-' _L and 1/ = 

o Jn' J v'> Jn' 2 

Pn pnpp Pnpp Pp 


2 2 

PppPn PnnPp PnpPnPp 

Pll — — P22 — ——, Pl2 — P21 — -To-, 

P P P 


P — PnnPpp PnpPpn- 


Note that pn + pi2 = Pn and P22 + P21 = Pp- 


( 12 ) 


(13) 

(14) 


If the equilibrium structure is axisymmetric about the rotation axis, the time dependence and 
(/)—dependence of the perturbations can be given by exp(icrt -f imcj)) with a being the oscillation frequency 
observed in an inertial frame, and m is an integer representing the azimuthal wavenumber. Introducing the 
vectors and defined as 





iujpn ’ 





iujpp ’ 


and 


iuj ’ 


(15) 


where lo = a + mfl is the oscillation frequency observed in a corotating frame of the star, the mass 
conservation equations are linearised to be 


P; + V • (pnC) = 0, 


(16) 
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and 

p;+V • ippe) = 0 . ( 17 ) 

Note that the mass conservation equation for the electron fluid becomes the same as that for the proton 
fluid because of the assumption of perfect charge neutrality, that is, and Pp/nip = p'^lmf^. The 

velocity equations (6) and (8) are linearized as 

^—F (D + —F = -V (/x; + vh ')+{e - e) x (v x uo), (is) 

Pn Pn P 

and 

—F (D +(— + —\f ie) = -V {'p' + {e - C) X (V X Uo), (19) 

Pp \pp rup) p^ 

where p = pp + peme/mp, and 

F (I) = —Lua^ + iuivo ■ + ioj^ ■ Vvg- (20) 

The Poisson equation is reduced to 

V2vI/' = 47rG(p; + Cp;). (21) 

Using a variant of Gibbs-Duhem relation, the pressure perturbation is given by 

P' = Pnp'n + PpP'j (22) 

where we have again ignored the entropy carried by the electron normal fluid. Note also that the superfluids 
carry no entropy. 

To obtain a relation between the densities and pp and the chemical potentials pn and p, we begin 
with writing the energy density e as 

e = e {pnj Pp )) (23) 

with which the chemical potentials are defined as 


pn {Pm Pp) — ^Pn )> Pp {Pni Pp) — i^^/^Pp)p^ ■ 

If the chemical potential of the electron is given by 


Pe{Pe) = (37r2;i3pe/m-eC3)2/3, 


we have, assuming pe = ppme/rup, 


(P'n \ ^ (Vll 

UV V^21 


where 


Vii 



Fi2 — 


We write the inverse of equation (26) as 


f dpn \ 
\dppj 


Pn 




( P'n \ ^ f 2ll 2 i2\ ( Pn\ 
\Pp ) \ 221 Q 22 ) \p'j' 


(24) 

(25) 

(26) 

(27) 

(28) 
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where 


(Qii 

Qi 2 \ 

(Vii 

Vl 2 \ 

VQ21 

Q22 J 

\V 21 

P22 J 


(29) 


To represent the rotationally deformed equipotential surfaces of a rotating star, we employ a coordinate 
system (a, 0, (j)), the relation of which to spherical polar coordinates (r, 9, (j)) is given by r = a[l + e(a, 0)], 
where e is proportional to and represents a small deviation of the equipotential surface from the 
corresponding spherical equipotential surface of the non-rotating star. We apply Chandrasekhar-Milne 
expansion (Chandrasekhar 1933a,b, Tassoul 1978) to the hydrostatic equations to determine the function 
e in the form e{a,9) = a{a) + P{a)P2{cos9) with P 2 being the Legendre function. See Lee (1993) for the 
definition of a{a) and f3{a). Since for uniformly rotating stars all the physical quantities in hydrostatic 
equilibrium are constant on the deformed equipotential surface, labeled by the coordinate a, we write the 
linearised basic equations using the coordinates {a,9,(j)) (Saio 1981, Lee 1993). In our formulation, the 
terms up to of order of in the perturbed velocity equations are retained so that the eigenfrequencies of 
the r-modes are correctly determined to the order of (see Yoshida & Lee 2000a). 


The perturbations in a uniformly rotating star are expanded in terms of spherical harmonic functions 
with different Vs for a given m (e.g., Lee & Saio 1986). For example, the vector is given by 


^ S]^{a)Yr{9,cl>)P‘ 


a , , 

/>|m| 

cn 

Si = 

/>|m| 


l>\m\ 




and the Euler perturbation of the gravitational potential, is given as 

vI/'= ^ vl/;(a)l^-(0,</,)e-‘, 

/>|m| 


(30) 

(31) 

(32) 

(33) 


where I = \m\ -I- 2{k — 1) and V = I + 1 for even modes, and I = \m\ + 2k — 1 and V = I — 1 for odd modes, 
and A: = 1, 2, 3, • • •. Substituting these expansions and the like into the linearised basic equations (16) 

~ (19) and (21), we obtain oscillation equations given as a set of simultaneous linear ordinary differential 
equations of the expansion coefficients (see Appendix), which is to be integrated in the superfluid core. 
The oscillation equations solved in the normal fluid envelope are the same as those given in Yoshida & Lee 
(2000a). 


To obtain a complete solution of an oscillation mode, solutions in the superfluid core and in the normal 
fluid envelope are matched at the interface between the two domains by imposing jump conditions given by 


L=C, U=ea, [p']t=0, ['l'']t=0, and [dd/Vda]! = 0, (34) 

where [/(a;)]l = lims^o{/(a^ + s) — f{x — s)}. The boundary conditions at the stellar center is the regularity 
condition of the perturbations g'/g, and 'i’'/g, where g = GMa/a^. The boundary conditions 

at the stellar surface are 6p = 0 and At' oc where Sp is the Lagrangian perturbation of the pressure. 

For numerical computation, oscillation equations of a finite dimension are obtained by disregarding the 
terms with I larger than Imax = |w| -I- 2kmax — 1 in the expansions of perturbations such as given by (30) to 
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(33). For the r-modes with V = \m\ calculated in this paper, we usually use kmax = 6 so that we can get 
reasonable convergence of the eigenfrequency and the eigenfunction. We solve the oscillation equations of 
a finite dimension as an eigenvalue problem of the oscillation frequency u using a Henyey type relaxation 
method (see, e.g., Unno et al 1989). 


3. Dissipations 


The stability of an oscillation mode of a star is determined by summing up all contributions from 
various damping and excitation mechanisms. If we consider the contributions from gravitational radiation 
reaction, viscous processes, and mutual friction in the superfluid core, the energy loss (or gain) rate dE/dt 
of a normal mode in a rotating neutron star may be given by 


dE 

dt 


= —aoj 


^ + \JlJ) - I d^x\^J2 + Cb IV • V'f 


1=2 


— 2n / d'^XPnBn 


dE 

dt 


GD 


dE 

dt 


PnPp 


GJ 


(w ■ w — ) 


dE\ 


dt J . 


dE 




dt 


2Cs 


fdE\ 


\dt J 


(35) 


MF 


where the asterisk (*) indicates the complex conjugate of the quantity, and the canonical energy E of 
oscillation observed in the corotating frame of the star is defined in the normal fluid envelope as (Friedman 
& Schutz 1978) 

VT' • 


FI = i y d^x ^pv' ■ v'* + —p'* — 
and in the superfluid core as (Mendell 1991b) 


47rG 


77’ / f f^.P / /* I ^ /T-s / /* 

E = - d X pv -v + —w w +2^ Vijp.pj - 


VT' • 


(36) 


(37) 


where v' = (/9„u(j + Ppv'^)/p, w' = v'p — v'^, p\ = p'^, and P 2 = p'p in the core. 

The terms {dE/dt)GD and {dE/dt)Gj on the right hand side of equation (35) denote the energy loss 
(or gain) rates due to gravitational radiation reaction associated with the mass multipole moment D'l^ and 
the mass current multipole moment where 


J' — 


and 


c{l + 1) 
Ni = 


DL = J d\p'r%"^*, 

J d^xr^ {pv' + p'vq) ■ (r x VY)™*), 
47rG (/ + !)(/+ 2) 


c2'+U(/-l)[(2l + l)!!]2’ 


(38) 

(39) 

(40) 


and c is the velocity of light (Thorne 1980, Lindblom et al 1998). 

The terms {dE/dt)s and {dE/dt)B are the energy dissipation rates due to shear and bulk viscosities, 
and Cs and (,b are the shear and the bulk viscosity coefficients, and crb is the traceless stress tensor for the 



perturbed velocity field (e.g., Landau & Lifshitz 1987). In this paper, we ignore the contribution from the 
bulk viscosity, which is important only for newly born hot neutron stars without superfiuids in the core. 
The shear viscosity coefficient we use in the superfluid core is 


Cs = 6 X 10 


18 


OX 


flO^KV 
~ ) 


g/cm s 


IQiSg/cm^^ 

(Cutler & Lindblom 1987, Sawyer 1989), and that in the normal fluid envelope is given by 


(41) 




(42) 


(Cutler & Lindblom 1987, Flowers & Itoh 1979). The stress tensor crC is evaluated by using v'^ in the 
superfiuid core (Lindblom & Mendell 2000). 


The term {dE/dt)MF is the energy loss rate due to mutual friction in the superfluid core, and the 
dimensionless coefficient is given by (Mendell 1991b) 


Bn = 0.011 X 


Pp_ ((/vV ( 


Pn\PpJ \PppJ \ 10i4gcm-3 


Pp 


1/6 


(43) 


Mutual friction is a dissipation mechanism inherent to a rotating system of superfiuids, and it is caused 
by scattering of normal fluid particles off the vortices in the superfiuids. Since we have assumed perfect 
charge neutrality between the electrons and protons, we consider scattering between the normal electrons 
and vortices of the neutron superfluid (Mendell 1991b, see also Alpar etal 1984). 


As is indicated by the first two terms on the right hand side of equation (35), if a normal mode has an 
oscillation frequency that satisfies —auj > 0, the oscillation energy E in the corotating frame, in the absence 
of other damping mechanisms, increases as a result of gravitational wave radiation, indicating instability 
of the mode (Friedman & Schutz 1978). It was Andersson (1998) and Friedman & Morsink (1998) who 
realized that the r-modes have oscillation frequencies that satisfy the instability condition. 


The damping (or growth) time-scale r of a normal mode may be given by 

1 _ 1 fdE\ _ 1 111 

T 2E \ dt J TGD tgj ts tmf ’ 


(44) 


where = 2E/{dE/dt)i. For the r-modes of I' = \m\, it is convenient to derive an extrapolation formula 
of the time-scale r given as a function of O and the interior temperature T (e.g., Lindblom et al 1998, 
Lindblom & Mendell 2000): 

r(fl,r) ^T^j\TrGp) V T ) ^ t1^p\ttGp) ’ 

where p = M//i), and only the dominant term in each of the dissipation processes with I = \m\ 4-1 = 3 
has been retained for the r-modes. The quantities Tq^, Tqj, t^, and r^p are assumed to be only weakly 
dependent on O and T. 


4. Numerical Results 

Following Lindblom & Mendell (2000), we employ a polytropic model of index A^ = 1 with the mass 
M = 1.4M0 and the radius R = 12.57km as a background model for modal analysis. The model is divided 
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into a superfluid core and a normal fluid envelope, the interface of which is set at p = ps = 2.8 x 10^‘^g/cm^. 
In the normal fluid envelope, we use the polytropic equation of state given hy p = Kp^ both for the 
equilibrium structure and for the oscillation equations. The core is assumed to be filled with neutron and 
proton superfluids and a normal fluid of electron, for which an equation of state, labeled A18+5v+UIX 
(Akmal, Pandharipande, and Ravenhall 1998), is used to give the energy density (23) and the relation (26) 
used for the oscillation equations. For the mass density coefficients Ppp, and in the core, we employ 
an empirical relation given by (see Lindblom & Mendell 2000) 

Ppjp « 0.031 + 8.8 X 10-i^p, (46) 


and a formula given by 


Pnp — PPnt 


(47) 


where rj is a parameter of order of ~ 0.04 that parametrizes the entrainment effects between the two 
superfluids (Borumand, Joynt, & Kluzniak 1996). 


We find that the r-modes oil' = m in the superfluid core are split into ordinary r-modes and superfluid 
r-modes, which we call r°-modes and r'^-modes. The toroidal components iTii of the r°- and r®-modes 
of /' = m = 2 are plotted versus a/R for the case oi rj = 0.04 and 17 = 17 /a/ GM/= 0.01 in Figure 
1, where the solid, dotted, dashed, and dash-dotted curves are used to indicate respectively the toroidal 
components iT”, iT^+ 2 , and iT ^+2 in the superfluid core, and the amplitude normalization is given 
by max(|zTm|) = 1. The figure shows that the two superfluids in the core flow together for the r°-modes, 
while they counter-move for the r^-modes. Note that the amplitudes of for the r^-mode are not 

necessarily negligibly small compared with those of iT^ in the core. This splitting of the r-modes in 
the superfluid core into two distinct families has been suggested by Andersson & Comer (2001). For the 
r°-mode, the radial dependence of the difference in the core is given in Figure 2, which shows 

that the difference is quite small compared with the normalization max(|zTm|) = 1. It is important to note 
that the I' = m r-modes having basically nodeless and dominant iT^ are the only r-modes we can find, as 
in the case of the r-modes in isentropic models (see Yoshida & Lee 2000a,b). 

In Table 1, the expansion coefficients kq and K 2 i and the scaled damping (or growth) timescales r? for 
the V = m = 2 r-modes are tabulated for the cases of ry = 0, 0.02, 0.04, and 0.06, where the coefficients kq 
and K 2 are defined in the expansion: 

w/f7 = ko(’ 7) + K2(?7)f^^ + 0(17"^). (48) 

Note that the exponent 7 employed to define in equation (45) is 7 = 2.5 for the r°-modes and 7 = 0.5 
for the r®-modes for the non-zero 77 ’s in the table. This is because the velocity difference w' = v'^ — v'^ 
approximately scales as w' oc 17^ for the r°-modes and as oc 17 for the r®-modes (see Lindblom & 
Mendell 2000). The coefficient kq for the r°-modes is numerically equal to 2m/[l'{l' + 1)], and K 2 is almost 
independent of the entrainment parameter rj. The value of kq is the same as the value found by Lindblom 
& Mendell (2000) and the value of K 2 differs only by 2.5%, suggesting that the r°-modes are the same 
modes found by Lindblom & Mendell (2000). (Because of different normalization conventions the value 
of K 2 reported by Lindblom & Mendell must be multiplied by a factor of 4/3 before comparing with our 
results.) The coefficients kq and K 2 for the r®-modes, on the other hand, appreciably depends on ry, and kq 
deviates from 2m/[l'{l' + 1)] as ry is increased from ry = 0. This kind of deviation of kq from 2m/[l'{l' + 1)] 
for the V = m r-modes has been found for relativistic neutron stars where the relativistic factor GM/c^R is 
regarded as a parameter (Yoshida 2001; Yoshida & Lee 2002a). Computing the r®-mode of F = to = 2 as a 
function of ry, we find that a linear formula given by 

Kg(ry) « 0.667-f 9.35ry 


(49) 
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gives a good fit to the r'*-mode frequency except at avoided crossings with inertial modes (see below). At 
r] = 0, the coefficients kq for the r°- and r'^-modes are both equal to 2m/[l'{l' + 1 )], which suggests that in 
the lowest order of the two r-modes are degenerate at 77 = 0. In Figure 3, the toroidal components iT^ 
and iTm +2 of the r°- and r®-modes of F = m = 2 at D = 0.01 are given versus a/R for the case of 77 = 0. 
The amplitudes of iTm +2 are much smaller than those of iT^ both for the r°- and r'^-modes. This figure 
also shows that \iT^ — 1 in the core for the r°-mode at 77 = 0 , for which we find that w' oc fl. 

Inertial modes in the superfluid core are also split into ordinary and superfluid inertial modes, which 
we call i°- and i®-modes. In Table 2, we have given kq for i°- and 7 ®-modes for m = 2 and 77 = 0, and see, 
e.g., Lockitch & Friedman (1999) and Yoshida & Lee (2000a) for the classification scheme employed here for 
inertial modes. For given m and Iq — |?ti|, we find at 77 = 0 pairs of i°- and t^-modes that have close values 
of kq, and the number of the pairs is equal to Iq — \m\. As an example, for the case of ttt = 2 and 77 = 0, 
the eigenfunctions iTii are shown for the i°- and z®-modes of kq = 0.5180 and kq = 0.5077 in Figure 4, for 
those of Kq = 0.4215 and kq = 0.4060 in Figure 5, and for those of kq = 1.1046 and kq = 1.1134 in Figure 6 . 
The inertial modes in Figure 4 belong to Iq — Ittt-I = 3 and those in Figures 5 and 6 to Iq — |w| = 5. Note 
that the r°- and r®-modes belong to Iq — |?ti| = 1 (e.g., Yoshida & Lee 2000a). It is generally observed for 
inertial modes with long radial wavelengths that the two superfluids co-move in the core for the 7 °-modes, 
and they counter-move for the 7 ®-modes. The coefficient kq for the 7 °-modes only weakly depends on the 
entrainment parameter 77 (see Figure 7). The coefficient kq for the z^-modes, on the other hand, increases 
approximately linearly as 77 is increased from 77 = 0 (see Figure 8 ). See Yoshida & Lee (2002b) for extended 
discussions on inertial modes in the superfluid core. 

The r®-modes (r°-modes) experience mode crossings with z°-modes (z®-modes) as the parameter 77 is 
increased from 77 = 0. For the case of ttz = 2 and fl = 0.01, Figure 7 illustrates an avoided crossing between 
the r^-mode of F = ttz and the z°-mode that tends to kq = 1.1046 as 77 ^ 0. In this figure, the dashed line 
is given by equation (49). For mode crossings between the r°-mode and z®-modes, on the other hand, it is 
quite difficult to numerically discern whether the mode crossings result in avoided crossing or degeneracy of 
the mode frequencies at the crossing point. Most prominent among such mode crossings of the r°-mode of 
V = m = 2 are those with the z®-modes that tend to kq = 0.5077 and kq = 0.4060 as 77 ^ 0. For the case of 
77 Z = 2 and fl = 0.01, the two mode crossings which occur at 77 « 0.0230 and 0.0484 are shown in Figure 8 , 
where the solid lines and the dashed line are for the z®-modes and the r°-mode, respectively. We note that 
the ratio w/fl « kq for the r°-mode is almost constant, equal to 2m/[V{l' + 1 )] = 0.6667 for V = m = 2, 
while the ratios w/fl for the z®-modes increase approximately linearly with increasing rj. Figure 9 shows the 
differences zT^ — iT^ versus a/R for the V = m = 2 r°-modes at 77 = 0.023 (panel a) and 0.0484 (panel b), 
and Figure 10 the eigenfunctions zT;/ for the z^-modes of ttz = 2 at 77 = 0.023 (panel a) and at 77 = 0.0484 
(panel b), where we have assumed Cl = 0.01 and the normalization max(|zTm|) = 1. The amplitudes of the 
differences at the mode crossing points are much larger than amplitudes of the difference off the crossing 
points (see, e.g., Figure 2). The resemblance between the z^-modes of ?o — |w| = 3 at 77 = 0 (Figure 4b) and 
at 77 = 0.023 (Figure 10a) and between the z®-modes oi Iq — \m\ = 5 at r] = 0 (Figure 5b) and 77 = 0.0484 
(Figure 10b) is obvious. 

7 ,From Table 1, we find that the growth timescales Tqj and Tqj^ for the r°-modes oil' = m = 2 are 
almost the same as those for the V = m = 2 r-modes of the N = \ polytropic model without superfluidity 
in the core (see, e.g., Yoshida & Lee 2000a). On the other hand, the growth timescales Tqj and Tqj^ for the 
r^-modes are much longer than those of the r°-modes, and the instability caused by the r®-modes is much 
weaker than the instability by the r°-modes. This is because the amount of gravitational wave radiation 
emitted from the r®-modes is much smaller than that from the r°-modes since (ppiT// + pniT//)/p ~ 0 and 
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p' ~ 0 for the r^-modes. It is interesting to note that Tqj of the r^-modes at p 0 is by several orders 
of magnitude longer than that at p = 0. Figure 11 shows, as a function of a/i?, 

for the r'^-modes at r/ = 0 (dashed line) and at p = 0.04 (solid line), where J'mmio) = dadJ^^/da. As 
found for the case of p = 0.04, the negative contributions to = R) the envelope almost completely 

cancel out the positive contributions in the core, which leads to extremely long growth timescales Tqj of 
the r®-modes at ry 0. 

Figure 12 illustrates the dependence of —T^p on the entrainment parameter rj for the V = m = 2 
r°-mode, and we find that —T^p has prominent and deep minimums at ry « 0.230 and 0.484., which is 
consistent with the result by Lindblom & Mendell (2000). Note that we have assumed ps = 2.8 x lO^^g/cm^ 
for the interface between the core and the envelope. These prominent minimums result from the mode 
crossings with the long radial wavelength i®-modes that tend to kq = 0.5077 and kq = 0.4060 as ry ^ 0 
(Figure 8; see also Andersson & Comer 20001). Mode crossings of the r°-mode with i®-modes that have 
much shorter radial wavelengths result in narrow and shallow minimums of —T^p, as found for ry >0.05. 
These minor dips were not found by Lindblom & Mendell (2000). We guess that the r-modes calculated by 
Lindblom & Mendell (2000) are less strongly affected by inertial modes associated with large Iq — \m\ having 
short radial wavelengths since they ignored terms higher than in their perturbative method. Comparing 
our results with those by Lindblom & Mendell (2000), we find that the damping timescale —T^p at a given 
ry is about by an order of magnitude shorter than that they obtained, and that our values of ry « 0.0230 and 
0.0484 for the local minimums of —T^p are slightly larger than their critical values 0.02294 and 0.04817. 
From Table 1 we also find that the value of /t 2 found here for the r°-modes is about by 2.5% larger than 
the value they found. We guess that these differences between the two calculations partly come from the 
differences in the way of evaluating the derivatives of the thermodynamical quantities that appear in the 
oscillation equations. For example, Lindblom & Mendell (2000) used p' = {dp/dp)pp' + {dp/d(})p(i', for 
which they assumed {dp/dp)p = {p/2p) from the polytropic relation p = Kp^ and employed a fitting 
formula to {dp/dj3)p obtained from Akmal et al’s equation of state (1998), where (i = p — In this 
paper, on the other hand, we used p' = (Qn + CQ2i)Mn + ( 2 i 2 + CQ 22 )p'p, and the coefficients Qij were all 
calculated by using equations (23) to (29) with Akmal et al’s equation of state (1998). 

The eigenfunction S(3{r, 9, </)) = dPo{r, 9, </)) + (4/3)5/32(?', 9, + 0(11'^) in Lindblom & Mendell (2000) 

may be given in terms of the eigenfunctions ^ and as 

6f3{r,0,(j)) = {Ma/M){R/a)ri~‘^'^Ay2^k{a)Yi^rni0,(l>), (50) 

k>l 

where r = a[l + e(a, 0)], and Aj/ 2 ,fc(a) = ki^) ~ 2^2 fc(®) Appendix for definition of the functions 

y 2 ^i/\a). Assuming 6Po = 0, we obtain 

d(32{r,9,^) « 0.75 x (M,/M)(i?/r)n-4 ^ Ayy2.fc(r)/i,™P,™(cos0)e*’"^ = ^ 5/32.fcP^^(cosd)e™^ (51) 

fc>l k>l 

where the mean radial distance a have been replaced by the radial distance r, and the factor fim is 
defined by the relation Y/^{9,<j)) = /imC)™(cos0)e*’"‘^, and Ik = \m\ + 2fc — 1. In Figure 13, we plot the 
functions Sf32,k{'>') for the r°-mode of /' = m = 2 at ry = 0.04, applying amplitude normalization given by 
y 2 ,k=i{R) = at the surface, where the solid, dotted, and dashed lines are for 5/32,i, 100 x 5/32,2, 

and 100 x 5/32,3, respectively. Since |5/32,i(r)| >> |5/32,fc(r’)| for k > 2, the 9 depenence of the function 
SP2{r,9,4>) is well represented by a single associated Legendre function P///pi{cos 9), and is not necessarily 
the same 9 dependence of the function 5/32 found by Lindblom & Mendell (2000). The amplitude of the 
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function 5(52s{r) is about by a factor of 3 larger than the amplitude of calculated by Lindblom & 
Mendell (2000), which is consistent with the result that T^p in this paper is about by an order of magnitude 
shorter than T^p they obtained. 


5. Discussion 


If we employ a set of the dependent variables defined as (see Lindblom & Mendell 2000) 


+ r = -, = C + P' = 

lUjp lUJ p 


(52) 


where p = pn + Pp, the continuity equations (16) and (17) and the velocity equations (18) and (19) are 
rewritten as 

Pn + Pp + ^- m = 0, (53) 


^_^+|.Vln^ 


P 


Pp 


Pn Pn Pp 


• V In ^ + V • = 0, 

P 


F (^) = - VI/+ ^ f V In 
F {O + ioj^^D iC) = -V/3', 

PnPp 

and the linearized Poisson equation remains the same: 

V2vl/' = 47rG {p'„+p'p), 


(54) 

(55) 

(56) 

(57) 


where D{^) = ^ x (V x Vq), and the terms of order of rrielmp have been ignored for simplicity. The 
term proportional to D {C") in equation (56) represents the drag force between the two superfluids. In a 
perturbative treatment of the r-modes, we may expand the eigenfunctions and eigenfrequencies in terms of 
fl as 


^ = ^o + ^2^^ + oin^), C = ^o+^2^^ + o{n^), (58) 

[/= + 0(0®), /3' = /3'02+/3^0‘‘ + 0(0®), (59) 

and 

w = kqO + K 2 O® + 0(0 ®), (T = sqO + S 2 O® + 0(0®), (60) 

where |^g| ~ \iTm\ ~ 0(1). Note that the quantities like V • ^ • Vp, Pp, and 4>' are of order of O^ in 

the lowest order for the r-modes. For the expansions given above, we have 

F(0 = F2 (Io, «0, So)0^ + F4(|o, ^2 , «0, «2, SO, 82 ) 0 ^ + 0(0®), (61) 


and we may not need to give the full expressions of F 2 and F 4 here. Equations (55) and (56) suggest 
the existence of two families of r-modes, that is, r-modes, which are governed by equation (55) and have 
Ko = 2m/[l'{l' + 1)], and r-modes, which are governed by equation (56) and have kq ^ 2m/\l'{l' -\- 1)]. If we 
assume ^ 2 (^ 0 ) = —VI /2 in equation (55), we obtain the lowest order r-mode solution given by iTm oc 
and Kg = 2m/[l'{l' + 1)], which leads to = 0 from equation (55), and to = 0 from equation (56). 
This suggests the existence of the r-mode solutions, for which kq = 2m/[l'{l' -I- 1)], and ^g yf 0 and U 2 yf 0, 
but = 0 and = 0. The solutions of this kind correspond to the r-modes Lindblom & Mendell (2000) 
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obtained in their perturbative treatment, and to the r°-modes found in this paper. On the other hand, 
because of the drag force, i.e., the term proportional to kq of the r-modes governed by equation 

(56) is not necessarily equal to 2m/[l'{l' + 1)]. For these solutions we may expect Ig ^ ^ ^2 ^ from 

equation (56), and 0 U 2 ^ 0 from equation (55). The solutions of this kind correspond to the 
r®-modes found in this paper (see also Andersson & Comer 2001). 

To show the importance of the drag force for the existence of the r-mode solutions with the scaling 
w' (X we calculate the r°- and r®-modes by ignoring the drag force terms proportional to D{^) on the 
right hand side of the velocity equations (18) and (19). The results are summarized in Table 3, in which the 
coefficients Kg, K 2 , and r°’s are tabulated for the r°- and r®-modes oi V = m = 2 for the case of 77 = 0.04. 
When we ignore the drag force terms, we can not find the r°-modes with the scaling of w' oc and 
the r-modes with the scaling w' (x Q are the only r-mode solutions we can find, and hence the exponent 
7 to define in the table is equal to 0.5 for both the r°- and r®-modes. In Figure 14, the toroidal 
components iT^ and iTm +2 of the two r-modes are given for rj = 0.04, where we have assumed = 0.01 
and the normalization maxdfTml) = 1. The amplitudes of iTm +2 are completely negligible compared to 
those of iTm- In the absence of the drag force. Kg’s for the two r-modes are both equal to 2m/[l'{l' + 1)], 
independent of rj, which means that the frequencies of the two r-modes are degenerate in the lowest order 
of O. These results suggest that the r-mode solutions with w' ex can be found only when the frequencies 
of the r-modes are not degenerate in the lowest order in fl. In this sense, it is reasonable to find to' oc O 
for the r°-modes at 7 = 0 since Kg’s of the r°- and r®-modes are equal to each other at 77 = 0, as shown by 
Table 1. 


6. Conclusion 

In this paper we have discussed the modal properties of the r-modes of neutron stars with the core 
filled with neutron and proton superfluids. We numerically find that the r-modes of rotating neutron 
stars with the superfluid core are split into ordinary r°-modes and superfluid r^-modes, and that the two 
superfluids in the core flow together for the r°-modes and they counter move for the r®-modes. These 
findings are consistent with earlier suggestions made analytically by Andersoon & Comer (2001). We also 
find that, although Kg of the r°-modes is numerically equal to 2m/[l'{l' + 1 )], almost independent of the 
entrainment parameter 77 , Kg of the r®-modes approximately linearly increases from 2m/[l'{V + 1)] as 77 is 
increased from 77 = 0. The r°-modes have the scaling of w' oc and are the same r-modes discussed by 
Lindblom & Mendell (2000), while the r^-modes have the scaling of w' oc fl. The instability caused by the 
r^-modes is found much weaker than the instability by the r°-mode and will be easily stabilized by various 
dissipation processes in the star. 

We have confirmed that, except in a few narrow regions of 77 around the prominent local minimums of 
T^p, the dissipation due to mutual friction in the superfluid core is ineffective to stabilize the instability 
by the r°-modes, as first shown by Lindblom & Mendell (2000). We have shown that these prominent local 
minimums of —T^p are caused by mode crossings between the r°-mode and the superfluid inertial z®-modes 
with long radial wavelengths comparable to those of the r°-mode (see also Andersson & Comer 2001). 
Since mutual friction is almost always very strong for the r®-modes and the instability by the r^-modes 
itself is quite weak, only the instability by the r°-modes will be of direct observational importance, e.g., as 
mechanisms that generate gravitational waves from normal modes of rotating neutron stars and/or cause 
the clustering of spin frequencies around 300Hz for neutron stars in LMXB’s, unless there exist other strong 
damping mechanisms for the r-modes (see Jones 2001a,b; Lindblom & Owen 2002; Haensel, Levenfish & 
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Yakovlev 2002). 

In this paper, we have made a brief report on a numerical result of inertial modes in the superfluid 
core. We find that inertial modes in superfluid neutron stars are also split into ordinary inertial i°-modes 
and superfluid inertial I'^-modes. It is generally observed for inertial modes with long radial wavelengths 
that the two superfluids co-move in the core for the i°-modes, and they counter-move for the i®-modes 
(see Yoshida & Lee 2002b for more complete discussions on inertial modes in the superfluid core). Inertial 
modes are expected to work as a mechanism limiting the amplitude growth of the r-modes (Morsink 2002). 

There are many interesting and challenging problems we have to deal with before we can conclude 
definitely about the neutron star oscillations. Confronting their cooling calculations of neutron stars with 
observational data for the surface temperature of several neutron stars, Kaminker, Haensel, & Yakovlev 
(2001), and Kaminker, Yakovlev, & Gnedin (2002) suggested that neutron superfluidity in the core of 
middle-aged neutron stars should be weak in the sense that the critical temperature Tc is less than 10®K. 
If this is the case, there exist no neutron superfluids in the core of temperatures T > 10®K. If neutrons in 
the core are in normal state, modal properties of low frequency oscillations propagating in the core will be 
different from those for the core filled with neutron and proton superfluids, since buoyant force in the core, 
produced by thermal and/or chemical stratification (Reisenegger & Goldreich 1992), comes into play. The 
normal neutron fluid core with stratification will support 5 -mode propagation, and the buoyant force in it 
will affect inertial modes in slowly rotating neutron stars. The nodeless r-modes of V = m, however, will 
remain almost the same even in the presence of buoyant force (see Yoshida & Lee 2000b), and damping 
due to mutual friction in the core will remain weak for the r°-modes, for which neutron and proton fluids 
co-move. Let us point out a possibility of using neutron star binaries as a probe to investigate existence or 
non-existence of neutron superfludity in the core, since tidally excited low frequency modes will be 5 -modes 
in the normal fluid core but inertial modes in the superfluid core, and the different tidal responses will result 
in different binary evolutions. The presence of a solid crust and a magnetic field in neutron star interior is 
another factor that makes complicated the problems of global oscillations of the stars. The solid crust is 
not a rigid body and supports its own normal modes (see, e.g., McDermott et al 1988, Lee & Strohmayer 
1996). For example, there exist torsional modes propagating in the crust, which can be resonantly coupled 
with the r-modes in the normal fluid core (Yoshida & Lee 2001; Levin & Ushomirsky 2001). Quite recently, 
using a local analysis of perturbations, Kinney & Mendell (2002) have suggested that no r-mode solution 
to the magnetohydrodynamic equations (e.g., Mendell 1998) exists in the superfluid core when both the 
neutron and proton vortices are pinned to the solid crust. This may suggest that a careful formulation 
of boundary conditions at the core-crust interfaces will be necessary to obtain reliable solutions to global 
oscillations of neutron stars, particularly when a magnetic field is essential for the oscillations. 


A. Oscillation Equations in the Superfluid Core 


For the oscillation equations in the superfluid core, we employ vectors y”, ^ 2 : I/ij y 2 J 
whose components are given by 


n _ on n _ 

VlM — ) i/ 2 .fc — 


4'; 


ga 


vi. = sr, 


p _ Bl 
y2,k — 




ga 


% 

y3,k = —, 

ga 


= (Al) 


y4,k = - , 

5 da 


where I = \m\ + 2{k — 1) for even modes and I = \m\ -I- 2A: — 1 for odd modes, and fc = 1, 2, 3, • • •. We also 
introduce vectors h", h^, and the components of which are 


K = Hr, = 


K = Hf, 


K = iTf,, 


(A2) 
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where V = I -\-1 for even modes and V = I — 1 for odd modes. Using these vectors, the perturbed continuity 
equations (16) and (17) are written as 


%i _ ( d\npn ^ dxsja) dxsjP) ^ 

da \ d In a da da 

+ (Aq + 3x3(/3)So) ^" + 3mx3(/3)(5o*^”) 


y^-—Qiiy^-—Qi2yl 

Pn Pn 


ga 


(Qll + C2 i2) 1/3 


(A3) 


dy\ 

‘ da 


dlnp. 


■P o dx^ia) dxsifd) , 

— o — (2-:- — (X -:-/if 


d In a da da 

(Ao + 3 x 3 (/ 3 )Bo) hP + ‘imx^{li)QoitP. 


yl-—Q2iy^-—Q22yl 


Pp 


Pp 


Pp 


{Q21 + C222) 2/3 


(A4) 


The radial components of equations (18) and (19) are reduced to 


dyo 

''~d^ 


ciw^—Uoy” + (1 - [7)2/2 + ciw^—U 02/1 

Pn Pn 

CIUJ^ ( —3/3^0 + mvE^ /i" - ciw2^3/3Boft.^ 

\ Pn J Pn 

ciw^ I ^ii3m/3(5o + vEiCq \ it" — ^^imfdQoitP, 

\ Pn ) Pn 


(A5) 


dyP 


da 


^ = ciw^ ( ^ j E^yP^ + (1 - [7)2/2 + ciU)'^—Eoyi 


P21 


CiW^ 


CiW^ 


Pp 

P22 ^ rUe 
Pp TTlp 


SPBo + mvC^Ei 


Pp 

hP - cico'^^SPBoh'" 
Pp 


^ + — ) 3m/3Qo + yt^EiCo 


Pp Tflp 

The perturbed Poisson equation (21) is reduced to 

dy-i 


itP — ciu;^^^3m/3Qo7t"- 
Pp 




(A6) 


(A7) 


— 47 rG'a^ (Qll + CQ21)^2/2 + 4:7rGa^ (Q12 + C222) 4 

+ I ^ [Aq — 47rGa^ (Qii + C22i + C(2i2 + QQ22)) 1 — 2 (al + / 3 Ao) Aq] 2/3 

+ / ^ [(—[7 + xo(ck)) 4 + Xo(/3)Ao — 6/3(Ao + Bo)] 2/4- (AS) 

Algebraic equations that determine the relations between the variables h, it, y^, and 2/2 are derived by 
making use of the horizontal components of the velocity equations ( 18 ) and ( 19 ), and they are 

- me - f„mFo{hP - le) - jm+iitp - le) = jy^ + /„j+(2/? - 2/1) + hi^ + (A 9 ) 

ClLU"^ ClUJ"^ 

LohP - mtP + fpmFoihP - h^) + fpM+iitP - le) = JyP - fpJ+{yP - y^) + hi^ + 622^, (AlO) 

CiUJ^ CiUJ^ 

-urn + - fnM^{hP - - fnmF,{ie - ft") = -Ky^ - fnK+{yP - y^), (All) 

-MohP + LiitP + fpM+{hP - h^) + fpmFiiitP - zt") = -KyP + fpK+{yP - y”), (A 12 ) 
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where 


fn = 

C = 1 + 


Pnp C 
Pn ^ 
We 

Wp’ 


fp — 


Pnp 1 

Pp C 


^ = 1 _|_ PPP 

Ulp Pp 


and 


/ Pn P12! Pn 

IPp {P22I Pp) + {rne/mp) 

The functions Xi(q^)j X 2 (a), X 3 (c«), and Xo(<a) are defined as 


fbii 

bi2\ 

\-( 

U 21 

^22 / 



(A13) 

(A14) 


(A15) 


, , da / \ t \ n 

XiW = (^ + a—, X2(a}=2a + a—, X3(a}=3a + a—, 

da da da 


(A16) 


and 


, , da d 

Xo(a) = +« 1 - 

da da 



(A17) 


The matrices Eg, Ei, Fq, Fi, I, J, J+, K, K~^, Lq, Zg , Li, Lf, Mq, Mq, Mi, and are defined as 


Ao = (1 + 2xi(a)) 1 + 2xi(/3)Ao, Ai = (1 + X 2 (a))l + X2(/3)Ao, (A18) 

T’o = ^Aoi[(l +2a + 2/3)1+ I 2 / 3 A 0 ], T’l = ^Aj-^ [(1 + 2a + 2/3)1 + 12/3Ai], (A19) 

/ = l-2xi(a)l-2xi(/3)Ao, (A20) 

J= J+-3/3Ao1(2Ao + Bo), (A21) 

J+ = mvAo^ [(1 + X2(a))l + x2(/3)Ao], (A22) 

K = k+ -3mf3ApQi, (A23) 

A'+ = i/A^ ^ [((1 + X2(a))l + x2(/3)Ai)Ci + 2((1 + X2(a) - X2(/3))l + 2x2{fd)Ai)Qi], (A24) 

Lq = (1 + 2a)l — w-Fq + /3Aq ^(2AoAo + QBq), (A25) 

Li = (1 + 2a)l - wF’i + f3Ap{2AiAi + 6 B 1 ), (A26) 

Mo = M+ -6mPApQi, (A27) 

M+ = ipAplil + 2a - 2/3)1 + 4f3Ai]QiAo + FiCi, (A28) 

Ml = M+ - QmPApQo, (A29) 

M+ = J^A-i[(l + 2a - 2/3)1 + 4/3Ao]QoAi + FqCo, (A30) 


where 1 denotes the unit matrix, and the matrices Aq, Ai, Bq, Bi, Cq, Ci, Qq, Qi, Aq, and Ai as well as 
the quantities U and ci are defined in Yoshida & Lee (2000a). 


The oscillation equations in the superfluid core are given as a set of linear ordinary differential equations 
for the variables with j = 1 to 4, which are obtained by eliminating the vectors h and it in equations 
(A3) to (A8) using the algebraic equations (A9) to (A12). 
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Table 1. 

Ko, K2- 

, and 

T°’s for the 

r°- and r'*- 

modes of V = 

m = 2. 

Modes 7 ] 

Kq 

K2 


'^GD (®) 



r° 

0 

0.6667 

0.408 

3.31 

3.83 X 10® 

-7.55 X 10'^ 

—oo 

T* ^ 


0.6667 

0.511 

2.12 X 10^ 

2.18 X 10® 

-2.99 X 10® 

—oo 

r° 

0.02 

0.6667 

0.408 

3.31 

3.83 X 10® 

-6.76 X 10^ 

-2.52 X 10® 

rs 


0.8532 

0.652 

4.76 X 10® 

1.21 X 10® 

-1.99 X 10® 

-2.94 X 10-1 


0.04 

0.6667 

0.408 

3.31 

3.83 X 10® 

-6.76 X 10'^ 

-1.70 X 10® 

r® 


1.0408 

0.787 

7.89 X 10® 

2.55 X 10® 

-1.14 X 10® 

-7.99 X 10“® 

r° 

0.06 

0.6667 

0.408 

3.31 

3.83 X 10® 

-6.76 X 10'^ 

-2.62 X 10® 

rS 


1.2296 

0.821 

3.64 X 10® 

1.03 X 10'^ 

-4.67 X lO"' 

-5.20 X 10“® 
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Table 2. kq for the i°- and i'*-modes for m = 2 and ij = 0. 


lo - \m\ 

i° 


2 

1.0978 

1.1190 


-0.5683 

-0.5109 

3 

1.3564 

1.3869 


0.5180 

0.5077 


-1.0293 

-0.9723 

4 

1.5191 

1.5526 


0.8639 

0.8612 


-0.2738 

-0.2417 


-1.2738 

-1.2238 

5 

1.6273 

1.6600 


1.1046 

1.1134 


0.4215 

0.4060 


-0.7028 

-0.6524 


-1.4340 

-1.3968 


Table 3. kq, K 2 , and r^’s for the r°- and r'^-modes oil' = m = 2 for rj = 0.04 without the drag force. 


Modes 

ACO 



'^GD 




0.6667 

0.408 

3.31 

3.83 X 102 

-7.49 X 10'*' 

-9.10 X 102 


0.6667 

0.529 

4.33 X 10"* 

2.63 X 10® 

-1.13 X 10® 

-7.30 X 10-2 
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Fig. 1.— Toroidal components iTm and iTm +2 of the r-modes of /' = m = 2 are plotted as a function of 
a/R for 0, = 0.01 and r] = 0.04, where the solid, dotted, dashed, and dash-dotted lines are used to indicate, 
respectively, iT”, iT”+ 2 > and iT ^+2 in the superfluid core. The amplitudes are normalized by the 
maximum value. Panels (a) and (b) are for the r°- and r'^-modes, respectively. 
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Fig. 2.— Difference iT]j^ — iT^ versus a/R for the V = m = 2 r°-mode for D = 0.01 and rj = 0.04, where the 
amplitude normalization rnaxdiTml) = 1 have been used. 

Fig. 3.— Same as Figure 1 but for the case of 77 = 0. 

Fig. 4.— Toroidal components iTm and iTm +2 of inertial modes of m = 2 and Iq — \m\ = 3 are plotted as a 
function of a/R for Cl = 0.01 and ?7 = 0 , where the solid, dotted, dashed, and dash-dotted lines are used to 
indicate, respectively, zT^, iT//, and in the superfluid core. The amplitudes are normalized by 

the maximum value. Panels (a) and (b) are for the z°-mode of kq = 0.5180 and the z^-mode of kq = 0.5077, 
respectively. 

Fig. 5.— Same as Figure 4 but for Iq — \m\ = 5. Panels (a) and (b) are for the z°-mode of kq = 0.4215 and 
the z®-mode of kq = 0.4060, respectively. 

Fig. 6 .— Same as Figure 4 but for Iq — \m\ = 5. Panels (a) and (b) are for the z°-mode of kq = 1.1046 and 
the z®-mode of Kq = 1.1134, respectively. 

Fig. 7.— Avoided crossing as a function of zy between the V = m r'^-mode and the z°-mode that tends to 
Kg = 1.1046 as ry ^ 0 , where m = 2 and Cl = 0.01 have been assumed. The dashed line indicates Kq given 
by equation (49) for the r'^-mode. 

Fig. 8 .— Mode crossings as a function of ry between the V = m r°-mode (dashed line) and the z®-modes 
(solid lines) that tend to kq = 0.5077 and kq = 0.4060 as ry ^ 0 , where m = 2 and Cl = 0.01 have been 
assumed. 

Fig. 9.— Differences iT// — iT// for the V = m = 2 r°-mode at D = 0.01 are plotted as a function of a/R 
for the cases of ry = 0.023 (panel a) and ry = 0.0484 (panel b). The amplitude normalization is given by 
maxdzTrnl) = 1 - 

Fig. 10.— Toroidal components zT^ and iTm +2 of inertial modes of rrz = 2 are plotted as a function of a/R 
for Cl = 0 . 01 , where the solid, dotted, dashed, and dash-dotted lines are used to indicate, respectively, iT^, 
iTm+ 2 , iT^, and iT ^+2 in the superffuid core. The amplitudes are normalized by the maximum value. Panel 
(a) shows the z^-mode of Iq — |zrz| = 3 at zy = 0.023, and panel (b) the z®-mode of Iq — |zrz| = 5 at zy = 0.0484 

Fig. 11.— NJ'mmio) versus a/R for the V = m = 2 r®-mode for the cases of ry = 0.04 (solid line) and zy = 0 
(dashed line), where N = a/|u;|A’ i/2A and J'mm{T) = Jq dadJmm/da- 

Fig. 12.- t'Iif versus zy for the /' = zrz = 2 r°-mode, where ps = 2.8 x lO^^g/cm^ has been used, and 

~'^MF i® given in second. 

Fig. 13.— Eigenfunction 5132 ,k as a function of r/R for the r°-mode of F = rrz = 2 for (7 = 0.01 and zy = 0.04, 
where amplitude normalizaation has been given by y 2 ,k=i(r = R) = fm+i,mCl^- See text for definition of 
the factor fim and the function y 2 ,k- The solid, dotted, and dashed lines are for 6 ^ 2 , 1 , 100 x 6 ^ 2 , 2 , and 
100 X 5^2,3, respectively. 
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Fig. 14.— Same as Figure 1 but for that we have ignored the drag force terms on the right hand side of the 
velocity equations (18) and (19) to calculate the r-modes. 
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